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BEAM  BREAKUP  INSTABILITIES  IN  LINEAR  ACCELERATORS: 
TRANSITION,  PHASE  MIXING  AND  NONLINEAR  FOCUSING 


The  control  of  beam  breakup  instabilities*  remains  a  major  concern  for 

7 

accelerator  builders.  Various  remedies  have  been  examined*":  lowering  the 

quality  factor  Q  of  the  deflecting  modes,  increasing  the  focal  strength, 

phase  mixing,  nonlinear  focusing,  and  detuning  of  the  cavities,  etc.  All 

of  these  may  readily  be  incorporated  in  the  continuum  models,* from  which 

asymptotic  formulas  have  been  obtained  for  certain  cases.  These  formulas 

have  frequently  been  used  to  compare  with  experimental  observations  and 
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with  other  more  elaborate  theoretical  works.  ’  At  the  moment,  they  are 

the  only  ones  available  to  describe  the  transition  of  BBU  from  the 

7  9 

cumulative  type  to  the  regenerative  type.  ’ 

In  spite  of  the  versatility  in  the  continuum  model,  a  direct  numerical 
integration  of  the  governing  partial  differential  equations  does  not  seem 
to  have  been  published,  although  the  evolution  of  BBU  has  been  calculated 
in  great  detail  by  other  means,  in  some  limiting  cases  .*'"*’ ***  Here,  we 
analyze  the  evolution  of  the  cumulative  BBU  by  a  direct  numerical 
integration,  intending  to  look  into  the  effects  of  phase  mixing  and 
nonlinear  focusing  which  defies  analytic  treatment.  This  work  was  also 
motivated  by  the  recent  observation^  that  the  asymptotic  growth  of  the 
cumulative  BBU  is  independent  of  the  linear  transverse  focal  strength.  It 
is  then  of  interest  to  determine  in  what  way  nonlinear  focusing  affects  the 
temporal  and  spatial  evolution  of  BBU. 

To  isolate  the  mechanism  of  phase  mixing  and  of  nonlinear  focusing,  we 
ignore  the  damping  associated  with  finite  Q.  Setting  Q  equal  to  infinity 
should  not  be  construed  as  this  damping  effect  being  unimportant.  It 
serves  the  additional  purpose  of  providing  a  more  stringent  test  for  the 
numerical  scheme.  For  without  a  finite  Q,  the  solution  grows  much  more 
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rapidly,  making  the  accuracy  in  the  numerical  integration  more  demanding. 

In  any  case,  the  effect  of  finite  0  enters  only  through  the  well-known 
decay  factor  exp  (-«ot/2Q)  in  the  amplitude  of  the  solutions,  where  is 
the  breakup  mode  frequency. 

Our  numerical  results  reveal  the  following.  When  >]  is  infinite,  the 
cumulative  BBU  evolves  asymptotically  (at  a  fixed  distance  from  the 
injector)  at  a  rate  which  is  independent  of  the  linear  transverse  focal 
strength.2  The  excellent  agreement  with  the  analytic  formula  adds  to  our 
confidence  in  our  numerical  calculation.  [In  a  separate  test,  the 
numerical  solution  agrees  with  the  analytical  results  to  within  20%  over 
nine  orders  of  magnitude  variations  in  the  beam's  transverse  displacement.] 
When  the  analysis  is  extended  to  include  a  finite  spread  in  the  betatron 
frequency,  we  find  that  this  finite  spread  has  virtually  no  effect  on  the 
above-mentioned  asymptotic  behavior.  While  neither  a  linear  focusing  force 
nor  a  finite  spread  in  the  betatron  frequency,  by  themselves,  can  provide 
long-term  stabilization  of  BBU,  a  sufficiently  strong  nonlinear  focusing 
force  may  limit  the  amplitude  of  the  beam's  transverse  displacement,  even 
in  the  limit  Q  -»  ®. 

For  a  continuous  coasting  beam,  the  evolution  of  the  cumulative  BBU  is 

1  3 

governed  by  the  following  normalized  equations  ’ 

(3/3T  +  3/3Z)2X  +  22X  =  A  (1) 

32A/3T2  +  A  «  2eX.  (2) 

Here  T  =  «  t(>0),  Z  =  «  z/c(>0),  2  =  k0c/M  where  ka  is  the  betatron 
o  o  pop 

wavenumber,  s  is  the  dimensionless  coupling  constant2  which  is  proportional 
to  the  beam  current  and  shunt  impedance  but  is  independent  of  the  focal 


2 


strength,  Q(Z,T)  is  a  measure  of  the  strength  of  the  deflecting  mode  which 

drives  the  BBU  and  X(Z,T)  is  the  normalized  transverse  displacement  of  the 

beam.  In  all  numerical  calculations,  we  assume  initial  rest  condition  and 

homogeneous  boundary  conditions;  only  one  non-trivial  boundary  condition  is 

imposed  at  Z  =  0:  X(0,T)  =  1  for  0  <  T  <  0.2  and  is  zero  otherwise.  Thus, 

X  is  normalized  to  the  initial  displacement,  imparted  upon  the  beam  at 

1  3-8 

Z  =  0  for  a  short  time  (similar  to  a  delta  function  excitation.  ’  ) 

Shown  in  Fig.l  is  the  spatial  distribution  of  X(Z,T)  at  four  different 
times,  obtained  from  a  fourth  order  numerical  scheme  in  space  and  time.  We 
set  e  =  0.000413  and  2  =  0.14.  The  dotted  curves  in  this  figure  are 
obtained  from  the  asymptotic  solution2: 

lX(Z’T)  I  “  B  3/4  exP  {  PCO  J  (3) 

where  x  -  e(T-Z)/ (S^Z)  >  0,  B  =  C/[Max(l,t/2.81) ]1/12, 


P(t) 


1.1573 


.1/6 


U  -  0.0556T  +  0.004x 


for  t  >  5 

for  t  <  5 


(4a, b) 


and  C  is  a  constant  to  be  determined  by  a  "one-point  fit"  with  the 
numerical  solution.  Equation  (3)  describes  the  transition  from  the  strong  / 


8  1 
focusing  solution  (p(t)  ->  1,  Eq.  (4b))  to  the  weak  focusing  solution  [Eq. 


(4a)]  as  T  increases.  This  transition  occurs  when  t  =  0.5.  The  excellent 
agreement,  over  a  wide  range  of  t,  between  the  numerical  solution  and  the 
analytic  theory  shown  in  Fig.  1  may  be  regarded  as  a  validation  for  both 


\  1 1'  i 

’  or-  ' 

V  v 


the  numerical  scheme  and  the  analytic  result  that  X  evolves  at  a  later  time 


as  if  the  linear  focusing  force  were  absent.  The  numerical  data  shows 


or 

s/ 
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that  the  spatial  variation  of  X  is  biased  toward  the  "slow  beam-cyclotron 
mode" . 2 

The  effect  of  phase  mixing  due  to  a  spread  in  the  betatron  frequency 
may  be  examined  by  pretending  that  the  electron  beam  consists  of  N  groups 
of  particles,  the  ith  group  bearing  its  individual  betatron  frequency  2. . 
Thus,  we  replace  Eq.  (1)  with  N  equations 

(3/3T  +  3/3Z)2X.  ♦  2i2Xi  »  A,  i  -  1,  ...  N  (5) 

and  the  quantity  X  in  the  right-hand  side  of  Eq.  (2)  is  replaced  by  <X>  = 

G^Xj  +  ^>2^2  +  •..■*■  G^XN  where  G^  represents  the  fraction  occupied  by  the 

114 

ith  group  (G^  +  G£  +  ...  +  G^  =  1) .  In  the  "spread  mass"  model,  ’  R.  = 
1/2  2 

2(i/N)  and  G^  =  6i(N  -  i)/N(w  -  1).  The  solution  <X>  is  shown  in  Fig. 

2  for  s  -  0.000413,  2  =  0.14,  N  =  8.  A  virtually  identical  figure  was 
obtained  when  we  increased  N  to  16.  From  this  figure,  we  see  that  a  finite 
spread  in  the  betatron  frequency  has  no  effect  on  the  long-time  growth  of 
the  cumulative  BBU.  This  is  due  to  the  fact  that  the  cumulative  BBU  grows 
asymptotically  at  a  rate  independent  of  the  linear  focusing  force,  and  a 
finite  spread  in  the  betatron  wavelength  cannot  alter  this  long-time 
behavior. 

A  nonlinear  octupole  focusing  force  may  be  handled  also  by  solving 

2 

Eqs.  (1)  and  (2),  except  that  the  2  X  term  in  Eq.  (1)  is  now  replaced  by 
3 

AX  in  a  simplified  model.  The  evolution  of  X  is  shown  in  Fig.  3  for  the 
case  A  =  0.0392.  Here,  we  see  that  the  amplitude  of  X  is  limited  at  large 
T,  in  marked  contrast  to  the  cases  of  linear  focusing  and  phase  mixing.  It 
is  possible  that  the  change  in  the  phase  relationship^  between  A  and  X  due 
to  nonlinearity  might  have  contributed  to  the  stabilization  shown  in  Fig. 

3,  but  this  issue  remains  to  be  examined  in  the  future. 


4 


In  summary,  an  accurate  numerical  scheme  was  devised  to  analyze  the 
classic  model  of  cumulative  BBU.  It  is  suggested  that  this  line  of  attack 
could  be  fruitfully  applied  to  a  great  variety  of  problems  in  a  rather 
straightforward  manner.  Interesting  areas  such  as  phase  mixing  and 
nonlinear  focusing  were  analyzed  in  this  initial  study,  and  some  unexpected 
features  reported. 

Helpful  discussions  with  S.  Zalesak  and  D.  Chernin  are  gratefully 
acknowledged.  Ve  also  thank  M.  Friedman,  C.  V.  Roberson  and  P.  Sprangle 
for  their  support.  This  work  was  supported  by  the  Department  of  Energy, 
Contract  No.  DE-AI05-86-ER13585  and  by  the  Office  of  Naval  Research. 
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